Selective filtering of a dithered image for the purpose of inverse dithering

ABSTRACT

A method for selecting a filter from among a plurality of ordered filters within an inverse dithering system is disclosed. In response to receipt of a number of filter parameters associated with the plurality of ordered filters and a selected portion of a dithered image, a determination is made of the region of support for the plurality of filters from the filter parameters. From these regions of support, a lowest ordered filter having a region of support that excludes an image edge in the selected portion of the dithered image is selected.

BACKGROUND OF THE INVENTION

1. Technical Field of the Invention

The present invention relates to inverse dithering, and more particularly, to a method for selectively filtering a dithered image for the purpose of inverse dithering the dithered image.

2. Description of Related Art

Dithering has long been a method to maintain the perceptual quality of continuous tone images in low bit depth video displays. The illusion of true continuous tone images can be produced on display devices incapable of providing a required number of amplitude levels necessary for achieving a full, continuous tone image by dithering. By dithering the images they may be displayed on systems having low bit-depth frame buffers while maintaining a reasonable image quality. The technique is also referred to as half-toning when the result produces a binary (black and white) image. Dithering enables a respectable display at a lower amplitude resolution by distributing quantization noise over small areas.

One process for improving the quality of a dithered image is referred to as inverse dithering. Inverse dithering strives to restore the original continuous tone nature of a predithered image. In inverse dithering, low-pass filters are used to reverse the effects of dithering such as graininess of the picture. However, use of low-pass filtering tends to blur the object edges of the image. Therefore, the inverse dithering process must ultimately achieve two conflicting goals. First, the effects of the original dithering must be smoothed, and at the same time the object edges or the high-frequency content of the original image must be maintained without undesirable blurring.

Inverse dithering may be discussed in a larger context of an image-rendering system. As illustrated in FIG. 1, a simplified image-rendering system consists of a dithering system 10 that receives as input, pixel values of an image I(x,y). The originally received image may be a monochromatic or a color image. The original image I(x,y) will have an amplitude resolution of q″ bits per pixel and the pixel values will be within the range of

I(x,y)∈{0,1,2 , . . . , 2^(q′)−1}.  (1)

The image I(x,y) is reduced in amplitude resolution to m bits per pixel by the dithering process. The dithered image I_(d)(x,y) will consist of pixels having values within the range of

 I_(d)(x,y)∈{0,1,2, . . . , 2^(m)−1}.  (2)

The dithered image I_(d)(x,y) is processed by an image buffer 15 which in many computer display systems would typically comprise a frame buffer. However, the image buffer 15 may comprise any number of other apparatus for holding or transporting an image such as a communication channel. The image buffer 15 actually drives the necessity for amplitude reduction of the original image I(x,y) due to the limited amount of memory in the buffer for storage or the narrow bandwidth of a communication channel.

The inverse dithering system 20 receives the dithered image I_(d)(x,y) and produces a reconstructed image I_(r)(x,Y) having a higher amplitude resolution of q″-bits per pixel. The amplitude resolution of the original and reconstructed images may be the same, but this is not required. The reconstructed image will have pixel values within the range

I_(r)(x,y)∈{0,1,2, . . . , 2^(q″)−1}.  (3)

A variety of methods have been developed for performing inverse dithering. The majority of these methods have been focused on the specialized case of inverse half-toning or, in other words, recovery of full gray scale from binary images. In one technique, a gray scale image is constructed from a binary image by means of a statistically generated lookup table. This method utilizes statistics generated during the original dithering process to perform the statistical analysis necessary to generate the lookup table. The statistical analysis requires complex calculations, and the storage of the lookup table greatly increases the memory requirements.

Another method utilizes cascading of adaptive run-length algorithms with statistical smoothing and impulse removal. A class of iterative, nonlinear algorithms have been proposed that are based on the theory of Projection Onto Convey Sets (POCS). This technique involves performing multiple iterations of an inverse dithering process until a final inverse dithered image is determined. This requires a great deal of processing due to the multiple iterations. Other methods have used information generated by a wavelet decomposition to perform the inverse dithering process.

There are also methods for reconstructing a dithered color image referred to as “color recovery”. The color recovery technique requires a renormalization process that increases the overhead requirements for processing the inverse dithering algorithm.

Some of the aforementioned techniques require prior knowledge of the dither mask used to dither the original image. Also, most of these methods involve highly complex calculations and require a great deal of memory in order to perform the inverse dithering processes. These requirements limit the ease and speed in which an inverse dithering process may be performed. Thus, there has arisen a need for a method of inverse dithering that is substantially less computationally complex and requires less memory resources than those of methods presently utilized.

OBJECTS OF THE INVENTION

Accordingly, it is an object of the present invention to provide a method and apparatus for inverse dithering an image which has been dithered according to any of a plurality of known dithering methods.

It is also an object of the present invention to provide an inverse dithered image using a plurality of ordered filters such that the filters process the dithered image according to a predetermined order.

It is yet another object of the present invention that one of the plurality of ordered filters may be selected such that the region of support of the selected filter does not include an image edge which is associated with a selected portion of a dithered image to be processed.

It is still further an object of the present invention to provide further advantages and features, which will become more apparent to those skilled in the art from the disclosure, including the preferred embodiment, which has been described below.

SUMMARY OF THE INVENTION

The following invention overcomes the foregoing and other problems with a system and method for selecting a filter from among a plurality of ordered filters within an inverse dithering system. The system includes a plurality of digital filters organized in a preselected order. The digital filters are either horizontally symmetric digital filters or a member of a horizontally asymmetric digital filter pair. The filters are ordered according to a predetermined set of rules such that filters having a highest filter index never comprise a subset of a lower indexed filter. The filters are also ordered such that asymmetric digital filter pairs always include adjacent filter indices. Finally, digital filters covering the same region of support have adjacent indices and are ordered an increasing cutoff frequency order.

A processor or other data processing unit is configured using hardware, software, firmware, etc. to process a number of parameters associated with the plurality of digital filters and data describing a selected portion of a dithered image. The filter parameters are used to determine regions of support for each of the plurality of ordered filters. The regions of support comprise areas that each of the ordered filters may process during a filtering operation.

A selection routine is implemented to sequentially cycle through the plurality of ordered filters such that a filter is initially selected from the plurality of ordered filters in accordance with the preselected order. The preselected order involves selection of a filter having a lowest index value among the plurality of filters. An edge detection process is performed to determine whether the selected filter includes an image edge within the region of support of the selected filter. A portion of an image contains an image edge if it does not have a constant color or gray scale level.

If an image edge is included within the region of support of the presently selected filter, the next highest ordered filter (i.e., the next sequential filter index) is selected for processing. If the selected filter region of support does not include an image edge, a determination of whether the presently selected filter is a member of a asymmetric filter pair is made. If so, the region of support of the second member of the asymmetric filter pair is processed to determine if it includes an image edge. If the second member includes an image edge, the presently selected filter is used for processing of the selected portion of the dithered image. If an image edge is not included within the region of support of the second member of the asymmetric filter pair, an indication is provided that no filtering is to be done on the selected portion of the dithered image.

If an image edge is not included within the region of support of the presently selected filter and the filter is not a member of an asymmetric filter pair, a determination of whether the filter is a member of a group of substantially identical support filters is made. If not, the presently selected filter is used for processing of the selected portion of the dithered image. For a filter belonging to a group of substantially identical support filters, a determination is made of the number of the group of substantially identical support filters. Next, a computation of an activity measure of the region of support for the group of identical support filters is calculated and used to select a member of the group of identical support filters to filter the selected portion of the dithered image. The activity measure indicates the magnitude of local variation in the dithered image. The selected filter is then used to process the selected portion of the dithered image.

BRIEF DESCRIPTION OF THE DRAWINGS

For a more complete understanding of the present invention, reference is made to the following detailed description taken in conjunction with the accompanying drawings wherein:

FIG. 1 is a functional block diagram of a simplified image rendering system;

FIG. 2 is a block diagram of the inverse dithering system of the present invention;

FIG. 2a is a block diagram of a processor implementation of the inverse dithering system of the present invention;

FIG. 3 is an illustration of a region of support of a low-pass filter;

FIG. 4a is an illustration of a region of support of a symmetrical filter;

FIG. 4b is an illustration of regions of support for a pair of asymmetrical filters;

FIG. 5 illustrates the ordering of a filter set according to the present invention;

FIG. 6 is a block diagram of the initialization system;

FIG. 7 is a block diagram of an alternative embodiment of the initialization system;

FIG. 8 is a block diagram of yet another alternative embodiment of the initialization subsystem;

FIG. 9 is a block diagram of the runtime subsystem;

FIG. 10 is an illustration of a windowed portion of an image;

FIG. 11 is an illustration of a flow chart for the filter selection process;

FIG. 12 is a block diagram of the filter and dequantization module;

FIG. 13 is a block diagram of an alternative embodiment of the filter and dequantization module;

FIG. 14 is an illustration of one embodiment of the up-multiplication module; and

FIG. 15 is a block diagram of a system for inverse dithering color images.

DETAILED DESCRIPTION OF THE INVENTION

Referring now to the Drawings, and more particularly to FIGS. 2 and 2A there is illustrated a block diagram of an inverse dithering system 20. The system 20 receives a dithered image I_(d)(x,y) in serial form, processes it and outputs a reconstructed image I_(r)(x,y). The inverse dithering system 20 is optimally implemented within a processor 21 utilizing software or firmware to perform the functionalities of the initialization subsystem 30 and runtime subsystem 35. The processor 21 receives the dithered image as input from video memory 22 and outputs the inverse dithered image to a display device 23. The inverse dithering system 20 works in conjunction with a filter set 25 including multiple filters 26 that enable the inverse dithering system to process the dithered image using an edge sensitive low-pass filtering approach. Each filter 26 of the low-pass filter set 25 include filter coefficients 28 which are provided to an initialization system 30 of the inverse dithering system 20.

The initialization subsystem 30 computes various parameters 32 required by the runtime subsystem 35 to generate the reconstructed image I_(r)(x,y). These parameters include the size of a window processed with each pixel of the dithered image, and, optionally, a correction lookup table utilized by the filtering process and, also optionally, a pattern match lookup table for use with a filter selection process. The runtime subsystem 35 utilizes the parameters provided by the initialization subsystem 30 to process the dithered image I_(d)(x,y) and generate the reconstructed image I_(r)(x,y).

The filter set 25 comprises a predefined set of digital low-pass, finite impulse response filters 26. Each filter 26 is defined by a number of filter coefficients and an impulse response, h_(i)(x,y), where (x,y) denotes the spatial coordinates relative to the pixel being processed and i comprises the index of the filter being used. Referring now to FIG. 3, there is illustrated the impulse response of a filter 26. The filter is defined as being finite-extent meaning that

h_(i)(x,y)=0 for x<−X_(i1) or x>X_(i2) or y<Y_(i1) or x>Y_(i2).  (4)

Equation (4) illustrates that the filter presented in FIG. 3 has an impulse response h_(i)(x,y) equal to zero for anything outside the region of support (R_(i)). The filter extents −X_(i1), X_(i2), −Y_(i1) and Y_(i2) represent the boundaries of the impulse response of the filter 26 on the X and Y axis.

Each point 42 within the region of support R_(i) 40 has associated with it a particular value representing the filter coefficient h_(i)(x₀,y₀) associated with that point (x₀,y₀). To maintain normalization, the filter coefficients h_(i)(x₀,y₀) (values) associated with each of these points 42 will sum to unity. $\begin{matrix} {{\sum\limits_{x}{\sum\limits_{y}{h_{i}\left( {x,y} \right)}}} = {1.0.}} & (5) \end{matrix}$

The region of support R_(i) for the filter is defined as

R_(i)={(x,y)|h_(i)(x′,y)≠0∀x′≦x and y′≦y}  (6)

and encloses a plurality of values h_(i)(x,y) denoting the filter coefficients within R_(i). The frequency response H_(i)(f_(x),f_(y)) of each filter 26 is given by $\begin{matrix} {{H_{i}\left( {f_{x},f_{y}} \right)} = {K{\sum\limits_{x}{\sum\limits_{y}{{h_{i}\left( {x,y} \right)}^{{- {j2\pi}}\quad {xf}_{x}}^{{- {j2\pi}}\quad {yf}_{y}}}}}}} & (7) \end{matrix}$

Given the frequency response the cutoff frequency ƒ_(c) of the filter 26 is defined as $\begin{matrix} {{\overset{\int\int}{\sqrt{f_{x}^{2} + f_{y}^{2}}} \leq {f_{c}{H_{i}\left( {f_{x},f_{y}} \right)}}} = {\frac{3}{4}\underset{{all}{({f_{x},f_{y}})}}{\int\int}{H_{i}\left( {f_{x},f_{y}} \right)}}} & (9) \end{matrix}$

The filters 26 within the filter set 25 are designed to fit into one of two categories. The filter 26 may be horizontally symmetric. In which case, the impulse response of the filter must have the following characteristics:

X_(i1)=X_(i2) and h_(i)(−x₀,y₀) for all (x⁰,y₀)  (10)

The region of support of a horizontally symmetric filter is represented in FIG. 4a where it can be seen that for every point on the left side of the Y axis there is a corresponding point on the right side of the Y axis.

Additionally, the filter 26 may comprise one member of a horizontally asymmetric filter pair in which case the impulse response of a filter pair h_(i)(x,y) and h_(j)(x,y) would have the characteristics:

 X_(i1)=X_(j2), X_(j1)=X_(i2) and h_(i)(−x₀,y₀)=h_(i)(−x₀,y₀) for all (x₀,y₀)  (11)

FIG. 4b illustrates a horizontally asymmetric filter pair wherein for every point on the left side of the Y axis for filter 50 there is a corresponding point on the right side of the Y axis in filter 55, and for every point on the right side of the Y axis for filter 50 there is a corresponding point on the left side of the Y axis in filter 55.

In addition to designing the filters 26 according to the above-identified constraints, the filters of the filter set 25 must be ordered according to a set of predetermined rules in order that the filter selection process may be more efficiently implemented within the runtime subsystem 35. The filters 26 are ordered according to the filter index (i). The filter index (i) is assumed to start at 1 and increase without interruption in increments of 1. The index 0 is reserved for the choice of no filtering. The filter index represents the order in which the filters will be considered during the filter selection process which will be more fully discussed in a moment.

The ordering of the filters 26 is subject to the following rules. First, in order for the filter selection process to make sense, filters 26 having lower index numbers must have an area of support R_(i) that extends in at least one direction further than the area of support R_(i) of filters having higher index numbers. In other words, the region of support of a filter 26 of index i cannot be a subset of the region of support of a filter having an index less than i. Second, asymmetric filter pairs as illustrated in FIG. 4b must have adjacent indices. This causes asymmetric filter pairs to be consecutively considered during the filter selection process. Finally, symmetric filters having the same region of support, which are referred to as “identical-support filters”, must also have adjacent indices. The identical-support filters will be ordered according to increasing cutoff frequencies. Thus, an identical-support filter with an filter index of 3 will have a lower cutoff frequency than an identical support filter having an index of 4.

Referring now to FIG. 5, there are illustrated the regions of support R_(i) for one particular instance of a filter set 25 consisting of six separate filters 26. The first rule for filter ordering is illustrated wherein it can be seen that the first filter having an index of 1 covers a larger region of support than any of the subsequent filters for indices i=2 through i=6. For each filter 26 the previous filters always covers some region not covered by the subsequent filter. The second rule is illustrated between filters 26 represented by indices i=3 and i=4. These filters 26 comprise an asymmetric filter pair having adjacent indices (i.e., 3 and 4). Finally, the filter pair having indices of 5 and 6 comprises and adjacent identical support filter pair wherein the filter having an index of 5 will have a lower cutoff frequency than the filter having an index of 6.

Referring now to the following example, there is illustrated one embodiment for a filter set 25 of the present invention, including five filters having various areas of support.

h₁(−4,0)=0.0156, h₁(−3,0)=0.1016, h₁(−2,0)=0.1406, h₁(−1,0)=0.1563, h₁(0,0)=0.1718,

h₁(1,0)=0.1563, h_(i)(2,0)=0.1406, h₁(3,0)=0.1016,

h₁(4,0)=0.0156,

h₁(x,y)=0.0 for x<−4, x>4, y<0, y>0.

h₂(−2,0)=0.0625, h₂(−1,0)=0.2500, h₂(0,0)=0.3750, h₂(1,0)=0.2500, h₂(2.0)=0.0625,

h₂(x,y)=0.0 for x<−2, x>2, y<0, and y>0.

h₃(−3,0)=0.0625, h₃(−2,0)=0.1250, h₃(−1,0)=0.3750, h₃(0,0)=0.4375,

h₃(x,y)=0.0 for x<−3, x>0, y<0, and y>0.

h₄(0,0)=0.04375, h₄(1,0)=0.3750, h₄(−2,0)=0.1250, h₄(3,0)=0.0625,

h₄(x,y)=0.0 for x<0, x>3, y<0, and y>0.

h₅(−1,0)=0.2500, h₅(0,)=0.5000, h₅(1,0)=0.2500,

h₅(x,y)=0.0 for x<−1, x<−1, x>1, y<0, and y>0.

Referring now to FIG. 6, there is illustrated a first embodiment of the initialization subsystem 30. The initialization subsystem 30 performs the processing that is necessary to generate the parameters and tables necessary for execution of the runtime subsystem 35 operations of the inverse dithering system 20. The initialization subsystem 30 routines are performed prior to the inverse dithering processes executed by the runtime subsystem 35. The initialization subsystem 30 receives filter coefficients 28 from the filter set 25. The filter coefficients 28 are input into an extraction of filter extents module 75 wherein the filter extents are extracted. The low pass filters are designed to have a finite impulse response as follows:

h_(i)(x,y)=0 for x<−X_(i1) or x>X_(i2) or Y>Y_(i1) or x>Y_(i2).  (12)

The values of X_(i1), X_(i2), Y_(i1), and Y_(i2) represent the extents of the filter 26 having impulse response h_(i)(x,y)=0 and can be derived from the filter coefficients 28.

The filter extents are input into the window size generator module 80 wherein the parameters for a sizing window for processing the dithered image are generated. The window size generator module 80 generates various parameters for a window to be used by the runtime subsystem 35. These parameters are defined as Δx, Δy₁ and ΔY₂. These parameters are dependent upon the extents of the filters 26 generated by the extraction of filter extents module 75. The window defined by the parameters must be large enough to cover the largest of the filters 26. The window parameters are given by the equations

Δx=max X_(i1)=max X_(i2),  (13)

Δy₁=max Y_(i1),  (14)

Δy₂=max Y_(i2).  (15)

The initialization subsystem 30 outputs the window parameters 85, filter extents 90 and filter coefficients 28 to the runtime subsystem 35 as the window parameter/lookup tables 32 illustrated previously with respect to FIG. 2.

Referring now to FIG. 7, there is illustrated an alternative embodiment for the initialization subsystem 30 wherein knowledge of the particular dither method utilized to originally process the dithered image is utilized to improve the inverse dithering process. For certain dither methods, edge detection in the inverse dithering module can be performed using pattern matching to achieve greater accuracy during the inverse dithering process. To enable the use of pattern matching to operate efficiently within the runtime subsystem 35, a precomputation of valid difference map patterns must be generated.

In the alternative embodiment illustrated in FIG. 7, the initialization subsystem 30 includes the filter extents extraction module 75 and window size generator 80 which operate exactly as previously discussed with respect to FIG. 6. Additionally, a pattern match lookup table generator module 95 is included for generating the pattern matching parameters. The pattern match lookup table generator module 95 utilizes a dither template 34 in conjunction with the extents of the filters 26 calculated by the filter extents extraction module 75 to generate possible dither patterns. From these generated dither patterns, valid difference maps of the region of support R_(i) of each filter 26 are derived. These are then forwarded to the runtime subsystem 35 in table form 96. Referring now to FIG. 8, there is illustrated yet another embodiment of the initialization subsystem 30, which utilizes implementation shortcuts for filtering. These shortcuts reduce the complexity of the system by using additional memory. One way in which the efficiency of the runtime subsystem 35 may be improved is to utilize a lookup table for correction terms used to implement various filtering modules which will be discussed more fully with respect to the runtime subsystem 35. This correction table may also be generated by the initialization subsystem 30. A correction term C_(ij) is generated for each filter and for each valid difference map pattern according to the equation $\begin{matrix} {C_{ij} = {G \cdot {\sum\limits_{{({u,v})} \in R_{1}}{{h_{i}\left( {u,v} \right)} \cdot {d_{ij}\left( {{x - u},{y - v}} \right)}}}}} & (16) \end{matrix}$

where i is the filter index, j is the index for the difference-map pattern, d_(ij) represents a valid pattern and G is the ideal gain represented by $\begin{matrix} {G = {\frac{2^{q} - 1}{2^{m} - 1}.}} & (17) \end{matrix}$

The table of correction terms is generated by the correction look-up table generator module 100. The module 100 outputs the terms in table form 101 for use by the runtime subsystem 35. The remainder of the modules within the embodiment illustrated in FIG. 8 operate as discussed previously.

Referring now to FIG. 9, there is illustrated a block diagram of the runtime subsystem 35. The runtime subsystem 35 includes a windowing module 105, filter selection module 110 and filtering and dequantization module 115. The runtime subsystem 35 receives as input the pixel values of a dithered image I_(d)(x,y) in serial form. This data is input to the windowing module 105. The windowing module 105 utilizes the dithered image and the window parameters 85 to generate a windowed portion of the dithered image for a presently processed pixel which is output to the filter selection module 110 and filtering dequantization module 115 over parallel interface lines 120.

The windowed portion comprises an area surrounding a presently processed pixel of the dithered image. Each pixel in the dithered image is processed with an associated windowed region until all pixels are processed. The filter selection module 110 selects a filter 26 most appropriate to filter the windowed portion of the image and provides the index of the filter to the filtering and dequantization module 115. The filtering and dequantization module 115 then utilizes the filter index to select the proper filter coefficients provided by the filter set 25 for processing of the windowed portion of the image and generates the pixel value of the reconstructed image I_(r)(x,y).

The windowing module 105 serves to highlight an area surrounding the particular pixel of the dithered image currently being processed by the runtime subsystem 35. Inverse dithering of an image is performed on a pixel-by-pixel basis. The pixels will be processed in the order in which they are generated or read from memory. Typically this means the pixels will be received in a serial raster fashion, i.e., the pixels are received from left-to-right and top-to-bottom. In theory, the order in which the pixels are received does not effect the result. However, the order in which the pixels are received does have some bearing on optimization of the process. The windowing process highlights a number of pixels surrounding the pixel currently being processed by the runtime subsystem 35. The size of the window highlighted around the currently processed pixel is given by the windowing parameters Δx, Δy₁, and Δy₂ provided by the initialization subsystem 30.

Referring now to FIG. 10, there is illustrated an example of the window 140 a generated around a pixel 145. The rectangular region comprising the window 140 a is defined by the four corners: (x−Δx, y−Δy₁), (x+Δx, y−Δy₁), (x+Δx,y+Δy₂), and (x−Δx, y+Δy₂). As mentioned previously with respect to the discussion of the initialization subsystem 30, the values of Δx, Δy₁ and Δy₂ are dependent upon the extents of the largest filter 26 within the filter set 25. Given that the pixel values in the region defined by the window 140 surrounding the processed pixel 145 are normally read into the windowing module 105 in a raster fashion, the windowing module 105 must contain enough memory to act as a buffer for storing the entire window portion 140 of the image and any pixels that would be read in between the window portions during a raster scan. In order to minimize the amount of memory required by the windowing module 105, Δy₁ and Δy₂ may be set to zero such that the windowed region 140 merely comprises a scan line including pixels on each side of the processed pixel 145.

In the case where the window region 140 b extends beyond the boundary of the image 135, as illustrated generally by the dashed rectangle 140 b, the values of the window region 140 b extending beyond the image 135 may merely be set to zero. However, it should be realized that there are other methods for dealing with a window region extending beyond the image 135 that are also consistent with the method and apparatus of the present invention.

The filter selection module 110 performs a filter selection based upon the values of the pixels within the window region 140. The filter selection module 110 operates according to a single constraint in that the filter selected by the module must only perform low-pass filtering on a region which at the original amplitude resolution comprises a constant color or constant gray scale level. In other words, the region of support of the selected low-pass filter 26 may not contain any edges of an image. Instead, the region of support may only include portions of the image having a constant color or gray scale level. Referring now to FIG. 11, there is illustrated a flow chart describing the filter selection process. The process involves the selection of the largest filter 20 from the filter set 25, which will fit the window region about the currently processed pixel that does not include an edge. The process begins at step 160 wherein the filter index i is set equal to 1. A filter index of i=1 correspond to the largest filter 26 within the filter set 25 as will be remembered from the previous discussion relating to ordering of the filters within the filter set. Next, at inquiry step 165, a determination is made if the filter index i equals the maximum filter index i. If so, this indicates that no filtering is to be performed on the presently processed pixel at step 170 and index zero is selected.

If inquiry 165 determines the filter is not the last filter, inquiry step 175 determines whether an edge of the image exits within the region of support (R_(i)) of the filter 26. If an image edge exist within the region of support, control passes to step 180 to increment the filter index i such that the next, smaller filter may be selected. Particular methods for determining whether an image edge exist within the region of support will be more fully discussed in a moment.

If an image edge is not detected in the region of support at inquiry step 175, inquiry step 185 determines whether the presently selected filter 26 is a member of an asymmetric filter pair. If so, inquiry step 190 determines whether an image edge exist within the region of support of the next filter within the filter set which comprises the asymmetric pair. If edges are absent in the regions of support of both asymmetric filters, the process returns at step 195 an index zero indicating no filtering of the pixel. Otherwise, the process returns at step 200 the current filter index.

If the current filter is not a member of an asymmetric filter pair, inquiry step 205 determines whether the filter is a member of an identical support group filter. If not, the current filter index is returned at step 210. If the current filter is a member of an identical support filter, the process will first determine the number of filters contained within the identical support filter group and set this equal to j at steps 215, 220 and 225. Once the number of filters within the identical support group is determined, the process calculates the activity measure A(R_(i)) of the region of support for the present filter according to the equation $\begin{matrix} {{A\left( R_{i} \right)} = {{\min \left\lbrack {\left( {\sum\limits_{{({u,v})} \in R_{1}}{{{I_{d}\left( {u,v} \right)} - {I_{d}\left( {x,y} \right)}}}} \right),\left( {E_{i} - {\sum\limits_{{({u,v})} \in R_{1}}{{{I_{d}\left( {u,v} \right)} - {I_{d}\left( {x,y} \right)}}}}} \right)} \right\rbrack}.}} & (18) \end{matrix}$

Using the activity measure, the (K+1)the filter in the identical support group will be selected at step 235 where K satisfies criteria $\begin{matrix} {\frac{j}{k} < \frac{2{A\left( R_{i} \right)}}{E_{i}} \leq {\frac{k + 1}{j}.}} & (19) \end{matrix}$

and can be rewritten as $\begin{matrix} {k = {{{ceiling}\left( \frac{2j\quad {A\left( R_{i} \right)}}{E_{i}} \right)} - 1.}} & (20) \end{matrix}$

Once k is calculated, the filter index (i+k) is returned at step 240. The activity measure indicates the magnitude of local variation in the dithered image.

The process of edge detection during the filter selection process described in FIG. 11 will now be more fully discussed. Edge detection is performed to insure that low-pass filtering is done only on regions of the image that are originally of a constant color or gray scale level. In a dithered image, constant color or gray scale regions are marked by the characteristic that the pixels values are always within one resolution level of each other. In some cases, there are recognizable patterns within the dithered region of an image. These particular characteristics lead is to edge detection methods that are unique and simple to implement.

One method of edge detection examines a map indicating differences between pixel values in a selected region (region of support) and the center pixel that is currently being processed. For the region to be declared edge free, the difference map must include either all zeros and ones are all zeros and negative ones. This indicates all pixels in the selected region are within one resolution level of each other. An example of this is illustrated in Table 1 below.

TABLE 1 Pixel Values Edge (center in bold) Difference Map Detection  5, 4, 5, 5, 4, 0, −1, 0, 0, −1 no edge 6, 7, 6, 6, 5 0, 1, 0, 0, −1 edge 5, 7, 6, 6, 6 0, 1, 0, 0, 0 no edge 5, 7, 4, 5, 4 1, 2, 0, 1, 0 edge

The region under consideration in line 1 includes pixels having resolution values 5,4,5,5,4 with the center pixel, having a value of 5, being the currently processed pixel. The difference map is generated at each position by subtracting the center pixel from each pixel position to generate the indicated difference map including 0, −1, 0, 0, −1. The difference map in the first line includes only zeros and negative ones. This indicates that no edge is included within the region. A similar indication is provided in line 3, wherein only zeros and ones are included within the difference map. Lines 2 and 4 include difference maps providing an indication that an edge exist within the region. In line 2, this is because both 1 and negative 1 are included with the zeros of the difference map. In line 4, this indication is provided by the 2 within the second position of the difference map.

The use of difference maps for edge detection methods is not a perfect process. In some cases, false negatives may be generated where the procedure fails to detect edges in the image. This leads to undesirable blurring of the restored image. For certain dithering methods, edge detection may be improved by observing that only a limited number of patterns can occur within the difference maps of regions having a constant color or gray level. Consider, for example, the dither patterns generated by an 8×1 one dimensional recursive tessellation array, such as [0,4,2,6,1,5,3,7]. The valid dither patterns for this array are illustrated in Table 2.

TABLE 2 Valid Dither Patterns 0, 0, 0, 0, 0, 0, 0, 0 0, 0, 0, 0, 0, 0, 0, 1 0, 0, 0, 1, 0, 0, 0, 1 0, 0, 0, 1, 0, 1, 0, 1 0, 1, 0, 1, 0, 1, 0, 1 0, 1, 0, 1, 0, 1, 1, 1 0, 1, 1, 1, 0, 1, 1, 1 0, 1, 1, 1, 1, 1, 1, 1

If the region of support is 5×1, the number of possible 5×1 difference maps be simply observed. For example,

All valid 5 by 1 difference maps are illustrated below in Table 3.

TABLE 3 Valid Difference Maps ±1, 0, 0, 0, 0 0, 0, 0, 0, ±1 0, ±1, 0, ±1, ±1 0, ±1, 0, 0, 0 ±1, 0, 0, 0, ±1 ±1, ±1, 0, ±1, 0 0, 0, 0, ±1, 0 0, ±1, 0, ±1, 0 ±1, ±1, 0, ±1, ±1,

Thus, once a difference map for a presently processed bit and region of support is determined, the difference map can be compared to difference maps within a table of valid difference maps indicating a constant color or gray scale at the filter selection module 110. If a corresponding difference map is found within the table then an edge does not occur within the presently processed region. If no corresponding valid difference map is detected an edge therefore must exist in the region of support. The difference map table is provided by the pattern match lookup table module 95 of the initialization subsystem 30.

If the phase by which the dither template was originally applied to the image is known, the number of valid difference maps at each specific pixel location may be further reduced. Using the more restrictive table of difference maps at each pixel location for pattern matching, it is possible to further reduce the number of false negatives arising in edge detection.

Referring now to FIG. 12, there is illustrated the filtering and dequantization module 115, which receives input bit values of the dithered image I_(d)(x,y), the filter coefficients 28 of the filter set 25 and the filter index of the selected filter 26 to generate an inverse dithered image I_(r)(x,y) having an increased bit depth. The received input pixel values I_(d)(x,y) have a resolution of m bits. These are up multiplied at the up-multiplied module 250 by gain G as defined in equation 17 to generate pixels of increased resolution with q-bit values. Therefore, the up-multiplied pixel values I_(s)(x,y) is given by

I_(x)(x,y)=G·I_(d)(x,y).  (21)

The up-multiplied pixel values I_(s)(x,y) are then processed by the filtering module 255 wherein the up-multiplied pixel values are multiplied and summed with the filter coefficients associated with the filter index provided by the filter selection module 210. The filtering is performed on the scaled pixels according to the equation $\begin{matrix} {{I_{r}\left( {x,y} \right)}{\sum\limits_{u}{\sum\limits_{v}{{h_{i}\left( {u,v} \right)}{{I_{s}\left( {{x - u},{y - v}} \right)}.}}}}} & (22) \end{matrix}$

Note that the first summation over u will contain X_(i1)+X_(i2)+1 terms whereas the second summation over V will contain Y_(i1)+Y_(i2)+1 terms. The filtering process may yield output values that are beyond the range of q-bits. When this happens, values above 2 ^(q)−1 are clamped to 2 ^(q)−1 while values below zero are clamped to 0.

Referring now to FIG. 13, there is illustrated an alternative embodiment of a more efficient implementation of the filter and the dequantization module 115. The filtering process described in the equation 22 above may be rewritten in the following manner. $\begin{matrix} \begin{matrix} {{I_{r}\left( {x,y} \right)} = \quad {\sum\limits_{u}{\sum\limits_{v}{{h_{i}\left( {u,v} \right)}{I_{s}\left( {{x - u},{y - v}} \right)}}}}} \\ {\quad {\sum\limits_{u}{\sum\limits_{v}{{h_{i}\left( {u,v} \right)} \cdot G \cdot {I_{d}\left( {{x - u},{y - v}} \right)}}}}} \\ {\quad {{\sum\limits_{u}{\sum\limits_{v}{{h_{i}\left( {u,v} \right)} \cdot G \cdot {I_{d}\left( {x,y} \right)}}}} + {\sum\limits_{u}{\sum\limits_{v}{{h_{i}\left( {u,v} \right)} \cdot G \cdot}}}}} \\ {\quad \left( {{I_{d}\left( {{x - u},{y - v}} \right)} - {I_{d}\left( {x,y} \right)}} \right)} \\ {\quad {{G \cdot {I_{d}\left( {x,y} \right)}} + {G \cdot {\sum\limits_{u}{\sum\limits_{v}{{h_{i}\left( {u,v} \right)} \cdot G \cdot \left( {{I_{d}\left( {{x - u},{y - v}} \right)} -} \right.}}}}}} \\ \left. \quad {I_{d}\left( {x,y} \right)} \right) \end{matrix} & (23) \end{matrix}$

The correction term C_(i)(x,y) may be defined as $\begin{matrix} {{C_{i}\left( {x,y} \right)} = {\sum\limits_{u}{\sum\limits_{v}{{h_{i}\left( {u,v} \right)}{G\left( {{I_{d}\left( {{x - u},{y - v}} \right)} - {I_{d}\left( {x,y} \right)}} \right)}}}}} & (24) \end{matrix}$

The terms (I(x−u,y−v)−I(x,y)) are the same as the entries of the difference maps described previously with respect to the discussions of edge protection methods. This greatly simplifies the computation of the correction term. A look-up table provided by the correction lookup table generator 100 of the initialization subsystem 35 may be used for its generation. Take, for example, the case wherein the 5×1 region of Table 3 there are only 18 possible cases. The look-up table for the 5×1 filter would only contain 18 entries. The dequantization and filtering module which increases the bit depth of the dithered image may now be described by the equation:

I_(r)(x,y)=G·I_(d)(x,y)+C_(i)(x,y).  (25)

Thus, the input dithered images are multiplied at up-multiplied module 250 and added together at adder 275 with the correction term generated by the difference map module 60, which generates a difference map of the processed bits, and the correction table 265 which compares the difference map from module 260 with a lookup table provided by the initialization subsystem 30. Note that we have only one multiplication as opposed to before.

Finally, as illustrated in FIG. 14, the reconstruction of a q-bit integer value from a m-bit integer value within the up-multiplied module 250 value can be approximated by a bit replication process. This eliminates the last remaining multiplication step. In the illustrated case a 3-bit integer value is scaled to an 8-bit integer value. This process is more fully described in copending U.S. application Ser. No. 09/015,935, which is incorporated herein by reference.

In the case of color images, each pixel is represented by a three dimensional vector wherein each vector corresponded to one of three channels: Red, green or blue, i.e,

I(x,y)=(I_(R)(x,y),I_(G)(x,y),I_(B)(x,y)).  (26)

or alternatively, may represent the luminance and the two chrominance channels as in

I(x,y)=(I_(y)(x,y),I_(u)(x,y),I_(v)(x,y))  (27)

Color images are processed by separating them into their individual components before inverse dithering each component. This is illustrated in FIG. 15, wherein a dithered image I(x,y) is provided to a channel separator 300 to separate the image into red, green and blue vectors 305. The red, green and blue vectors 305 are inverse dithered at 310 in the manner described previously. The inverse dithered images are then input to a combining module 315 wherein the color image is reconstructed from the inverse dithered vectors.

Although preferred embodiments of the method and apparatus of the present invention have been illustrated in the accompanying Drawings and described in the foregoing Detailed Description, it is understood that the invention is not limited to the embodiments disclosed, but is capable of numerous rearrangements, modifications, and substitutions without departing from the spirit of the invention as set forth and defined by the following claims. 

What is claimed is:
 1. A method for selecting from a plurality of ordered filters to filter dithered images, comprising the steps of: determining regions of support for a plurality of ordered filters from filter parameters for the plurality of ordered filters; selecting a lowest ordered filter having a region of support that excludes an image edge within a selected portion of the dithered image; and using the selected lowest ordered filter to filter the selected portion of the dithered image.
 2. The method of claim 1, wherein the step of selecting further includes the steps of: selecting a first lowest ordered filter of the plurality of ordered filters; determining whether the selected first filter includes the image edge within the region of support of the selected first filter; identifying the selected first filter as a selected lowest ordered filter only if the image edge is excluded from the region of support of the selected first filter; selecting a second lowest ordered filter of the plurality of ordered filter if the image edge is included in the region of support of the selected first filter; and identifying the selected second filter as the selected lowest ordered filter only if the image edge is excluded from the region of support of the selected second filter.
 3. The method of claim 2, wherein the first and the second lowest ordered filters are in a preselected order.
 4. The method of claim 1 further comprising the steps of: determining if a selected first filter is a member of an asymmetric filter pair only if the image edge is excluded from the region of support of the selected first filter; determining whether the image edge exists in a other member of the asymmetric filter pair if the selected first filter is a member of the asymmetric filter pair; identifying the selected first filter as the selected lowest ordered filter only if the image edge is included in the region of support of the other member of the asymmetric filter pair.
 5. The method of claim 4 further comprising the step of: indicating no filtering of the selected portion of the dithered image only if the image edge is excluded from the region of support of the other member of the asymmetric filter pair.
 6. The method of claim 1 further comprising the steps of: determining if a selected first filter is a member of a group of substantially identical filters; identifying the selected first filter as the selected lowest ordered filter only if the selected first filter is not the member of the group of identical support filters; determining which member of the group of identical support filters best filters the selected portion of the dithered image if the presently selected filter is the member of the group of identical support filters; and identifying the member which best filters the selected portion of the dithered image as the selected lowest ordered filter.
 7. The method of claim 6, wherein the step of determining which member of the group of identical support filters best filter the selected portion of the dithered image further includes the steps of: determining a number of members of the group of identical support filters; identifying an activity measure of a region of support of the group of identical support filters; and selecting the member of the group of substantially identical filters which best filter the selected portion of the dithered image based upon the activity measure.
 8. An article of manufacture for selecting between a plurality of ordered filters to filter a dithered image, comprising: computer readable storage medium; and computer programming stored on said storage medium; wherein said computer programming is configured to be readable from said computer readable storage medium by a processor and thereby causing said processor to operate so as to: determine a region of support for a plurality of ordered filters from filter parameters for the plurality of ordered filters; and select a lowest ordered filter having a region of support that excludes an image edge within a selected portion of the dithered image; and use the selected lowest ordered filter to filter the selected portion of the dithered image.
 9. The article of manufacture of claim 8, wherein the computer programming is further configured to cause said processor to: select a first lowest ordered filter of the plurality ordered filters; determine whether the selected first filter includes the image edge within the region of support of the selected first filter; identifying the selected first filter as the selected lowest ordered filter only if the image edge is excluded from the region of support of the selected first filter; select a second lowest ordered filter if the image edge is included in the region of support of the selected first filter; and identifying the selected second filter as the selected lowest ordered filter only if the image edge is excluded from the region of support of the selected second filter.
 10. The article of manufacture of claim 9, wherein the first and the second lowest ordered filters are in a preselected order.
 11. The article of manufacture of claim 8, wherein the computer program is further configured to cause said processor to: determine if a selected first filter is a member of an asymmetric filter pair only if the image edge is excluded from the region of support of the selected first filter; determine whether the image edge exists in an other member of the asymmetric filter pair if the selected first filter is the member of the asymmetric filter pair; identify the selected first filter as the selected lowest ordered filter only if the image edge is included in the region of support of the other member of the asymmetric filter pair.
 12. The article of manufacture of claim 11, wherein the computer program is further configured to: indicate no filtering of the selected portion of the dithered image only if the image edge is excluded from the region of support of the other member of the asymmetric filter pair.
 13. The article of manufacture of claim 8, wherein the computer programming is further configured to cause said computer to: determine if a selected first filter is a member of a group of substantially identical filters; identifying the selected first filter as the selected lowest ordered filter only if the selected first filter is not the member of the group of identical support group filters; determine which member of the group of identical support group filters best filters the selected portion of the dithered image if the presently selected filter is the member of the group of identical support filters; and identify the member which best filters the selected portion of the dithered image as the selected lowest ordered filter.
 14. The article of manufacture of claim 13, wherein the computer programming is further configured to cause said computer to: determine a number of members of the group of identical support filters; identify an activity measure of a region of support of the group of identical support filters; and select a member of the group of substantially identical filters which best filter the selected portion of the dithered image based upon the activity measure.
 15. A system for selecting a filter, comprising: a plurality of digital filters organized in a preselected order; and a processor configured to select a lowest ordered filter including a region of support that excludes an image edge within a selected portion of a dithered image.
 16. The system of claim 15, wherein the processor is further configured to: select a first lowest ordered filter; determine whether the selected first filter includes the image edge within the region of support of the selected first filter; identify the first selected filter as the selected lowest ordered filter only if the image edge is excluded from the region of support of the first selected filter; select a second lowest ordered filter if the image edge is included within the region of support of the selected filter; and identify the first selected filter as the selected lowest ordered filter only if the image edge is excluded from the region of support of the first selected filter.
 17. The system of claim 16, wherein the first and the second lowest ordered filters are in a preselected order.
 18. The system of claim 15, wherein the process is further configured: determine if a selected first filter is a member of an asymmetric filter pair only if the image edge is excluded from the region of support of the selected first filter; determine whether the image edge exists in an other member of the asymmetric filter pair if the selected first filter is the member of the asymmetric filter pair; identify the selected first filter as the selected lowest ordered filter only if the image edge is included in the region of support of the other member of the asymmetric filter pair.
 19. The system of claim 18, wherein the processor is further configured to: indicate no filtering of the selected portion of the dithered image only if the image edge is excluded from the region of support of the other member of the asymmetric filter pair.
 20. The system of claim 15, wherein the processor is further configured to: determine if a selected first filter is a member of a group of substantially identical filters; identifying the selected first filter as the selected lowest ordered filter only if the selected first filter is not the member of the group of identical support group filters; determine which member of the group of identical support group filters best filters the selected portion of the dithered image if the presently selected filter is the member of the group of identical support filters; and identify the member which best filters the selected portion of the dithered image as the selected lowest ordered filter.
 21. The system of claim 20 wherein the processor is further configured to: determine a number of members of the group of identical support filters; identify an activity measure of a region of support of the group of identical support filters; and select a member of the group of substantially identical filters which best filter the selected portion of the dithered image based upon the activity measure.
 22. The system of claim 15, wherein each of the digital filters of the plurality of digital filters comprises at least one of a horizontally symmetric digital filter and a member of a horizontal asymmetric digital filter pair.
 23. The system of claim 15, wherein the plurality of digital filters are ordered such that asymmetric filter pairs always have adjacent filter indices.
 24. The system of claim 15, wherein the plurality of digital filters are ordered such that filters having the same region of support have adjacent filter indices and are ordered in increasing cut-off frequency order. 